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ABSTRACT. We study interior regularity of solutions of a generalized stationary Stokes 
problem in the plane. The main, elliptic part of the problem is given in the form div(A(Du)), 
where D is the symmetric part of the gradient. The model case is A(Du) = (k+ Du|)''~ 2 Du. 
We show optimal BMO and Campanato estimates for A(Du). Some corollaries for the 
generalized stationary Navier-Stokes system and for its evolutionary variant are also men- 
tioned. 

1. Introduction 

Let Q. C R 2 be a domain. In this article we study properties of the local weak solution 
u € W l ' lf '(Q.) and n € Lf* (£1) of the generalized Stokes problem 1 

-divA(Du) + Vrc = -divG in ft, 

divu = in Ll 

for given G : Q. — > M 2 * 2 . Here u stands for the velocity of a fluid and n for its pres- 
sure. We do not need boundary condition, since our results are local. The model case is 
A(Q) = v(k+ |Q|) p ~ 2 Q corresponding to power law fluids with v>0, K" > 0, I < p <°° 
and Q symmetric. But we also allow more general growth conditions, which include for 
example Carreau type fluids A(Q) = H*>Q + v(k+ |Q|) ;, ~ 2 Q with fi^ > (see Subsec- 
tion 2.2). In this article we are interested in the qualitative properties of A(Du) and n in 
terms of G. The divergence form of the right-hand side is only for convenience of the 
formulation of the result, since every f can be written as — divG with G symmetric, see 
Remark 3.12. 

System (1.1) originates in fluid mechanics. It is a simplified stationary variant of the 
system 

(1.2) u r -divA(Du) + [Vu]u + V7r=-divG, divu = 0, 

where u stands for a velocity of a fluid and n for its pressure. The extra stress tensor 
A determines properties of the fluid and must be given by a constitutive law. If A(Q) = 
2vQ with constant viscosity v > 0, then (1.2) is the famous Navier-Stokes system, which 
describes the flow of a Newtonian fluids. In the case of Non-Newtonian fluids however, 
the viscosity is not constant but may depend non-linearly on Du. The power law fluids and 
the Carreau type fluids are such examples, which are widely used among engineers. For a 
more detailed discussion on the connection with mathematical modeling see e.g. [27, 28]. 
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'We denote by Du the symmetric part of the gradient of u, i.e. Du = (Vu + (Vu) r )/2. 
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The existence theory for such type of fluids was initiated by Ladyzhenskaya [24, 25] and 
Lions [26]. 

The main result of the article are the following Campanato type estimates for the local 
weak solutions of (1.1). 

Theorem 1.1. There is an a > such that for all £ [0, Ct) a C > exists such that for 
every ball B with 2B C Q. 

+R-P-f\A(Du) - (A(Du)) 2B |<foj 

2B 

In particular, G £ £? i2+ P(2B) implies A(Du), K £ _Sf l ' Z+ P(B). 

The spaces l2+ $(B) are the Campanato spaces, see Subsection 2.1. Our main theo- 
rem in particular includes the BMO-case (bounded mean oscillation), since BMO = jSf 1 ' 2 . 
Theorem 1.1 is a consequence of the refined BMO(B-estimates of Theorem 3.9, which also 
include the case VMO (vanishing mean oscillation). The upper bound a is given by the 
maximal (local) regularity of the homogeneous generalized Stokes system. Our estimates 
hold up to this regularity exponent. Due to the Campanato characterization of Holder 
spaces C 0,a our results can also be expressed in terms of Holder spaces. 

Our result is an extension of the results in [12] to the context of Non-Newtonian fluids. 
In [12] we studied the problem 

(1.3) -div(A(Vu)) = -divG 

for Q. C W, n > 1, with similar growth conditions on A but Du replaced by the full gra- 
dient Vu. The model case is A(Vu) = |Vu| p_2 Vu with 1 < p < °°. In [12] we proved 
Theorem 1.1 for weak solutions of (1.3). The case p > 2, has been studied first in [9]. 

Theorem 1.1 is the limit case of the nonlinear Calderon-Zygmund theory, which was 
initiated by [17, 18]. Iwaniec proved that G £ U with s > p' implies A(Vu) £ U, where 
p' = -^y. See also [9] for related works. In the context of fluids the corresponding result 
was obtained in [11]. However, the limiting regularity of (1.1) for G = restricts the 
transfer of integrability to the range s £ [p 1 , z^np'] for n > 3 and s £ [p' ,°°) for n = 2. 

The reduced regularity for (1.1) with G = is also the reason, why we can only treat 
the planar case n = 2 in this paper. The crucial ingredient for Theorem 1.1 are the decay 
estimates for the homogeneous case G = in terms of the gradients. In this paper we are 
able to prove such decay estimates in the planar case n = 2, see Theorem 3.8. If such 
estimates can be proven for « > 3, then Theorem 1.1 would directly generalize to this 
situation. Unfortunately, this is an open problem, even in the absence of the pressure. 

Theorem 1.1 can be used to improve the known regularity results for the stationary 
problem with convective term [Vu]u, see Section 5, and for the instationary problem (1.2), 
see Section 5. The first C l a -regularity for planar flows were obtained in the series of the 
articles [20, 21, 22] under various boundary conditions under the restriction K > 0. See 
also [30, 1]. The stationary degenerate case K > was treated in [32] for 1 < p < 2. To our 
knowledge the only result for n > 2 is the one obtained in [8] with K > and 1 < p < 2 and 
small data and zero boundary values. Because of the zero boundary values (combined with 
the small data), we are not able to use this result for the higher regularity of the case G = 0. 

Note that our result is optimal with respect to the regularity of G. All other planar 
results mentioned above need much stronger assumptions on the regularity of G. This is 
one of the advantages of the non-linear Calderon-Zygmund theory. This is the basis for our 
improved results in Section 4 and Section 5 for the system including the convective term. 
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It is based on the fact, that the convective term can be written as div(u ® u) using divu = 
and therefore can be treated as a force term divG. 

2. Notation, basic definitions and auxiliary results 

2.1. Notation. By B we will always denote a ball in R 2 . We write B <s £2 if the closure of 
B is contained in D.. Let \B\ denote the volume of B. Vector valued mappings are denoted 
by bold letters, e.g. u, while single valued functions with regular letters, e.g. 77. For 
f £ Lj (R 2 ) we define component- wise 

(f) B = U{x)dx := / t{x)dx, and M^f = / |f- (f) B \dx, (mH)(x) = supM^f. 

J \B\ Jb J B3x 

B B 

The space BMO of function of bounded mean oscillations is defined via the seminorm (for 
£2 open) 



Ibmo(o) : = SU P T l f - (?)B\dx = supM^f; 



saying that f G BMO(£2), whenever its seminorm is bounded. Therefore f 6 BMO(R 2 ) if 
and only if M#f e L°°(R 2 ). We say that a function f £ BMO(£2) belongs to the subspace 
VMO/ oc (n) if lim £ ^o + sup Bc£2 \g\ <£ M B f = 0. We need also the following refinements of 
BMO, see [31]. For a non-decreasing function co : (0,°°) — > (0,°°) we define 

Mlj=-^-J\f-(f)B\dx, 

B 

where R is the radius of B. We define the seminorm 

ll f llBMO m( n) : = su P M L,B f - 

sen 

The choice fi)(r) = 1 gives the usual BMO seminorm, while co(r) = r a with < a < 1 
induces the Campanato space «if 1-2+a (£2). Its seminorm we denote |u| jfi,2+fir a y 

By kB, with k > 0, we denote the ball with the same center and k times the radius. For 
functions f,g on Q. we define (f,g) := f n f(x)g(x)dx. Similarly also for mappings to R", 
n > 1. We write / ~ g if and only if there exist constants cq,ci > 0, such that 

cof<g<cif, 

where we always indicate on what the constants may depend. Furthermore, we use c, C (no 
index) as generic constants, i. e. their values may change line to line but does not depend 
on the important quantities. 

We say that a function p : [0,°o) — > [0,°°) is almost increasing if there is c > such that 
for all < s < t the inequality p (s) <cp(t) is valid. We say that p is almost decreasing if 
there is c > such that for all < s < t the inequality p(s) > cp(t) is valid. We say that p 
is almost monotone if it is almost increasing or almost decreasing. 

For a mapping u : D. -» R 2 we define Du = (Vu + (Vu) r )/2, Wu = (Vu - (Vu) r )/2 
and ([Vu]u) ; - = Li=i u fc<%U/- In the parts of the article dealing with evolutionary problems 
we will assume that u : £2 x (0 7 T) — >• R 2 . In this case all operators V, D, W and div are 
understood only with respect to the variable i£fl. 

For P, Q e R" with n > 1 we define P Q = £" =1 P ; Q;. The symbol R 2 ^ denotes the set 
of symmetric 2x2 matrices. For a set M C R" we denote %m as the characteristic function 
of the set M, i.e. %{ x ) = 1 if x £ M otherwise it is equal to zero. We write R-° = [0, +°°). 
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2.2. N-functions and the extra stress tensor A. The following definitions and results are 
standard in the context of N-function (see e. g. [29]). 

Definition 2.1. A real function <p : R-° — > R-° is said to be an N-function if it satisfies 
the following conditions: There exists the derivative <p' of (p. This derivative is right con- 
tinuous, non-decreasing and satisfies <p'(0) = and (p'(t) > for t > 0. Especially, (p is 
convex. 

Definition 2.2. We say that the N-function <p satisfies the ^-condition, if there exists 
c\ > such that for all t > it holds (p(2t) < c\ <p(f). By A2(<p) we denote the smallest 
constant c\. For a family <P of N-functions we define A2(<£ > ) := sup (f>e<1> A2((p). 

Let (p be an N-function. We state some of its basic properties. Since <p(f) < <p{2t) the 
A2-condition is equivalent to <p(2t ) ~ <p{t). The complementary function <p* is given by 

(p*(u) := sup {ut - <p(t)) . 

r>0 

It satisfies (<p*)'(f ) = (<p') _1 (t), where (<p')~ l is the right-continuous inverse of <p'. More- 
over, (<p*)* = (p. 

For all 8 > there exists c§ (only depending on A2(<p*)) such that for all t ,u > 

(2.1) tu < 8(p(t)+c S (p*(u). 
This inequality is called Young's inequality. For all t > 

(2.2) ^(0<9(O<^(O, 9 (^)<^)< 9 (^). 
Therefore, uniformly in f > 

(2.3) <p(r)~<p'(r)f, <p* ((?'(?)) ~ <p(t), 

where the constants only depend on A2((p,(p*). 

For an N-function (p with A2(<p) < °°, we denote by LP and W the classical Orlicz 
and Sobolev-Orlicz spaces, i.e. u € LP if and only if f (p(\u\)dx < °° and u G W 1<p if and 
only if u, Vu e Lf . By Wo' v (£l) we denote the closure of Cq (Q.) in W 1,p (£2). 

Throughout the paper we will assume that (p satisfies the following assumption. 

Assumption 2.3. Let (p be an N-function with (p e C 2 ((0, +°°)) (1C 1 ([0,+°°)) such that 
(p" is almost monotone on (0, +°°) and 

<p'{t)~t<p"{t) 

uniformly in t > 0. The constants hidden in ~ are called the characteristics of (p. 

We remark that if <p satisfies Assumption 2.3 below, then A2({<p,<p*}) < °° will be 
automatically satisfied, where A2({<p,<p*}) depends only on the characteristics of <p, see 
for example [2] for a proof. Most steps in our proof do not require that q>" is almost 
monotone. It is only needed in Theorem 3.7 for the derivation of the decay estimates of 
Theorem 3.8. 

Let us now state the assumptions on A. 

Assumption 2.4. Let <p hold Assumption 2.3. The vector field A : M 2x2 — > M 2x2 , A 6 
c 0,l( R 2x2 \ {0}) nC (M 2x2 ) satisfies the non-standard (p-growth condition, i. e. there are 
c,C > such that for all P,Q e R^ym with p 

(A(P) - A(Q)) • (P-Q) > c(p"(\P\ + |Q|) |P-Q| 2 , 
|A(P)-A(Q)| <C<p"(|P| + |Q|)|P-Q| 
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holds. We also require that A(D) is symmetric for all D £ IRgym and A(0) = 0. 

Let us provide a few typical examples. If <p satisfies Assumption 2.3. Then both 
A(Q) := <p'(|Q|)|§r and A(Q) := <p'(|Q sym |)jP^ satisfy Assumption 2.4. See [10] for 
a proof of this result. In this case, (1.1) is just the Euler-Lagrange equation of the local 
W^-minimizer of the energy ,X(w) := / Q <p(|Dw|)c/x + (G, Vw). Here wif is the sub- 
space of functions weff 1 '' 1 with divw = 0. The pressure acts as a Lagrange multiplier. 
This includes in particular the case of power law and Carreau type fluids: 

(a) Power law fluids with 1 < p < °°, K > and v > 

A(Q) = v(jc+ |Q|) P_2 Q and a>(t) = [' v(k + s) p - 2 sds 

Jo 

or 

A(Q) = v(jc 2 + |Q| 2 ) £ ^Q and <p(t) = f ' v(k 2 + s 2 )^ sds. 

Jo 

(b) Carreau type fluids with 1 < p < oo, K", ;Uoo > and v > 

A(Q)=M} + v(fc+|Q|) p ~ 2 Q and <p(t) = f ^s + v{K + s) p - 2 sds. 

Jo 

(c) For 1 < p < oo, jUoo > 0, and v > 

Q f 

A(Q) = /iooQ + varcsinh(|Q|)— T and (p(t) = / fXooS + V as:csmh(s) ds . 

|Q| Jo 

We introduce the family of shifted N-functions {<p a } a >o by <Pa(0A := (p'(a + t )/(a + t). 
If (jff satisfies Assumption 2.3 then <p"(f) ~ (p"(a +t) uniformly in a,t > 0. Moreover, 
(^fl)*})«>0 < °° depending only on the characteristics of <p. Let use define V : 

m>2x2 _v T»2x2 



(2-5) V(Q) = v /(p'(|Q|)|Q||g|. 

In the special case of A(Q) := <p'(IQI) jjjj ^ e quantity V(Q) is characterized by 

|V(Q„- A (Q,.Q a„a 

The connection between A, V, and the shifted N-functions is best reflected in the following 
lemma, which is a summary of Lemmas 3, 21, and 26 of [10]. 

Lemma 2.5. Let A. satisfy Assumption 2.4. Then 

(A(P)-A(Q)).(P-Q)~|V(P)-V(Q)| 2 

~<P|Q|(|P-Q|) 
-(«P*)|A(Q)|(I A ( P )- A (Q)D 

uniformly in P,Q £ M 2x2 . Moreover, 

A(Q)-QHV(Q)| 2 ~<p(|Q|), 

and 

|A(P)-A(Q)|~( 9 | Q| )'(|P-Q|), 

uniformly in P,Q £ K 2x2 . 



(1 



L. DIENING, P. KAPLICKY, AND S. SCHWARZACHER 



As a further consequence of Assumption 2.3 there exists 1 < p < q < °° and K\ > 
such that 

(2.6) <p(st) < Ki max{s p ,s q }q)(t) 

for all s,t > 0. The exponents p and q are called the lower and upper index of <p, respec- 
tively. We say that q> is of type T(p 1 q,K\) if it satisfies (2.6), where we allow 1 < p < q < °° 
in this definition. 

The following two lemmas show an important invariance in terms of shifts. 

Lemma 2.6 (Lemma 22, [13]). Let (p hold Assumption 2.3. Then (<p\p\)*(t) ~ (<P*)|a(p)|(0 
holds uniformly in t > and P G K 2x2 . The implicit constants depend on p, q and K\ only. 

We define 

(2.7) ~p := min{p,2} and q := max{<7,2}. 

Lemma 2.7. Let (p be of type T(p,q,Ki) and P G R 2x2 , then there is a K depending on 
K\,p,q such that <P|p| is of type T(p,q,K) and (<P\p\)* and (<P*)\a(p)\ are of type T(q' \~p ,K). 

The proof can be found in [12, Lemma 2.3]. Finally, we can deduce from Lemma 2.7 the 
following versions of Young's inequality. For all 8 G (0, 1] and all t,s > it holds 

ts<K^8 l -^(p(t) + 8(p*(s), 

(2.8) 

ts< 8(p{t) + K p - 1 8 1 -" <p*{s). 

3. Proof of the main theorem 

Let u, % be the local weak solution of (1.1), in the sense that u G W d if (H), n G if (£2), 
and 

(3.1) V<^ G W^(Q.) : (A(Du),D|> - (7r,div<^) = (G,D§), 

where we used that A(Du) and G are symmetric. To omit the pressure, we will use diver- 
gence free test function, i.e. 

(3.2) V<? G Wt;l(£l) : <A(Du),D§> = (G,DQ- 

The method of the proof of Theorem 1 . 1 was developed in [ 1 2] for elliptic systems with the 
main part depending on full gradient of solutions. It is based on a reverse Holder inequality, 
an approximation by the problem with zero right hand side and a decay estimate for this 
approximation. These three properties are discussed in the subsequent subsections. Note 
that the restriction to the planar case and <p" almost monotone is only needed for the decay 
estimate of Subsection 3.3. The first two subsections are valid independently of these extra 
assumptions. 

3.1. Reverse Holder inequality. In this section we show the reverse Holder estimate for 
solutions of (1.1). To prove the result we need a Sobolev-Poincare inequality in the Orlicz 
setting from [10, Lemma 7]. 

Theorem 3.1 (Sobolev-Poincare). Let y/ be an N-function such that y and y* satisfy the 
^-condition. Then there exists < 9 < 1 and c > 0, which depend only on A2({ \jf, V^*}) 
such that is almost convex 2 and the following holds. For every ballB C IR" with radius R 



It is proportional to a convex function. 
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and every v <G W lv (B) holds 



i 



(3.3) -f V ( ]y ^ }Bl \dx<cU V e (\Vy\)dx 

B B 

Remark 3.2. It is not possible to replace full gradient on the right hand side with the 
symmetric one only. Consider \ = (xx, — JCi) on the unit ball. 

We also need the following version of the Korn's inequality for Orlicz spaces, which is 
a minor modification of the one in [14, Theorem 6.13]. See [5] for sharp conditions for 
Korn's inequality on Orlicz spaces. 

Lemma 3.3. Let B CR" be a ball. Let \j/ be an N-function such that \j/ and satisfy the 
A2-condition (for example let \j/ satisfy Assumption 2.3). Then for all v G W^'^^B) with 
(Wv) b = the inequality 



V(|Vv|)flfx<C / y/(|Dv|)afx 

JB 

holds. The constant C > depends only on ^{{ V, W*}) < °°- 
Proof. From [14, Theorem 6.13] we know that 

(3.4) [ y(\Vv-(Vv) B \)dx<C [ \ir(\Dv-{Dv) B \)dx. 



Using (Wv)b = we have Vv = (Vv — (Vv)s) + (Dv)s. Thus, by triangle inequality 
and (3.4) we get 

\ir(\Vv\)dx<c [ \if(\Dv-{Dv} B \)dx + c f \jf(\{Dv) B \)dx, 

JB JB 

where we also used A2( v) < 00 ■ Now, the claim follows by triangle inequality and Jensen's 
inequality. □ 

As in [12] we need a reverse Holder estimate for the oscillation of the gradients. Addi- 
tional difficulties arise due to the symmetric gradient and the hidden pressure (so that the 
test functions must be divergence free). 

Lemma 3.4. Let u be a local weak solution of (1.1) and B be a ball satisfying 2B C fi. 
There exists 9 € (0,1) and c > only depending on the characteristics of <p, such that for 
all P, Go S IRjy^, 

i 



j |V(Du) - V{P)\ 2 dx <cij | V(Du) - V(P) \ 2e dx 

B 2B 

+ c-f(<P*)\ M p)\(\G-Go\)dx 



2B 

holds. The constant c > depends only on the characteristics of <p G T{p 1 q 1 K) and the 
constants in Assumption 2.4. 

Proof. Let 77 G Cq{2B) with % B < r] < Xib/2 aR d |Vt]| <c/R, where R is the radius of B. 
We define y/ = 77^(11 — z), where z is a linear function such that (u — z)?b = 0, Dz = P, and 
Wz = (Wu)2fl- We cannot use y/ as test function in the pressure free formulation (3.2), 
since its divergence does not vanish. Therefore, we correct \j/ by help of the Bogovskii 
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operator Bog from [3], In particular, w = Bog(divi//) is a special solution of the auxiliary 
problem 

divw = divy/ in jB 
w = md{\B). 

We extend w by zero outside of \B. It has been shown in [14, Theorem 6.6] that Vw can 
estimated by divi// in any suitable Orlicz spaces. In our case we use the following estimate 
in terms of <P| P | . 

j- 9| P |(|Dw|)<ix < C-j^ (p\ P \(\dw\j/\)dx. 

2B 2B 

The constant C > depends only on the characteristics of (p. 
Using divu = 0, we have 

divi// = (u - z) + r^div(u - z) = qrf~ l V7] (u - z) - ?7%P. 

This implies 

(3.5) ^(p| P |(|Dw|)<fcc < cj(p\ ¥ \ \^-Y) dx + C-l(p\ P \(\trP\)dx. 

2B 2B 2B 

We define £, := y — w = TJ^(u — z) — w, then div^ = 0, which ensures that ^ is a valid 
test function for (1.1). We get 

(3.6) 

(A(Du) - A(P), 77«(Du - P)) = (G - G , rj«(Du - P)> + (G - G , (u - z) ® sym V(ij*)) 

- (A(Du) - A(P), (U - Z) ®sym V(H ¥ )> 

- (G - G ,Dw) + (A(Du) -A(P),Dw). 

The symbol ® S ym denotes the symmetric part of <g>, i.e.(f <E>symg)ij '■= (fig; + f/g/)/2 for 
f,g £ R 2 . We divide (3.6) by \2B\ and estimate the two sides. Concerning the left hand 
side we find by Lemma 2.5 

\2Bf l (A(Du) - A(P),7f(Du -P)) ~ j jf | V(Du) - V(P) \ 2 dx =: (I). 

2B 

We estimate the right hand side of (3.6) by Young's inequality (2.8) for tpipi with 8 G (0, 1) 
using also (<P|p|)* ~ (<P*)|a(p)| ( see Lemma 2.6). 

(/) <c 5 ^(9*)|A(p)|(|G - G \)dx + 8j 77^(p|p| (|Du - ¥\)dx 

2B 2B 

+ csj(P\e\ y R Z j dx + c 5 j- 9 |P | (|Dw[ )dx 

2B 2B 

+ 5^ T ,9-i»'( 9 *) |A( p )I (|A(Du)-A(P)|) < fa 

2B 

=:(II) + (III) + (IV) + (V) + (VI). 
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Now we use Lemma 2.5 to estimate (///) + (VI) < 8 c(I), so these terms can be absorbed. 
Moreover, by (3.5) 



(IV) + (V)<c(IV)+c-f(p l p l (\tiP\)dx. 



IB 

Since P is constant, trP = divz and divu = 0, we can estimate 
(3.7) 

•^P[ P |(|trP|)dx= fy(<P|p|) e (|div(u-z)|VxJ < f ^(<P| P |) e (|Du - Dz)\)dx 

2B IB IB 

It remains to estimate (IV). We use Sobolev-Poincare inequality of Theorem 3.1 with \y = 
<P| P | such that (<P|p|) e is almost convex and 

(IV) =cjq> n dx < c U <P|||(|Vu - Vz\)dx\ * 

2B 2B 

with G (0,1). The constants and 9 are independent of |P|, since the A2({(p a } a>0 ) is 
bounded in terms of the characteristics of <p. 

As (W(u — z))2b = we find by Korn's inequality (Lemma 3.3) with y = <p^ (almost 
convex) and Dz = P that 



(IV) <c(j-tf v] (\T)u.-T)i\)dx 



IB 

The above estimates and Lemma 2.5 show that 



(IV) + (V) < c ^9 |p| (|Du-Dz|)^V <c^\V(Du)-V(P)\ 2e dx^j 



IB IB 

The lemma is proved. □ 

Lemma 3.4 allows to obtain exactly as in [12] the next corollary, compare with [12, 
Corollary 3.5]. 

Corollary 3.5. Let the assumptions of Lemma 3.4 be satisfied. Then for all P G Ksym 
j | V(Du) - Y(P)\ 2 dx < c (p*)|A(p)| (j |A(Du) - A(P)| dx 

B 2B 

+ c (<P*)\A(P)\(\\G\\bMO(2B)) 

The constants only depend on the characteristics of <p and the constants in Assumption 2.4. 

3.2. Approximation property. Let u be a local weak solution of (1.1) and B be a ball 
satisfying 2B C i2. We consider a solution h, p of the homogeneous problem 

-divA(Dh) + Vp = in H, 
(3.8) divu = in i2, 

h = u on d£l. 

The next lemma estimates the natural distance between u and its approximation h. 
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Lemma 3.6. For every 8 > there exists eg > 1 such that 

j |V(Du) - V(Dh)| 2 £ /x < 8 (<p*)|(A(Dn)>2s| (/ l A ( D «) - (A(Du)) 2B |^ 

B IB 

+ c 5(9*)KA(Du)> 2s |(l|G|| BMO(2B) ) 

holds. The constants depend only on the characteristics of (p and the constants in Assump- 
tion 2.4. 

Proof. The estimate is obtained by testing the difference of the equations for u and h by 
u h. The proof is exactly as in [12, Lemma 4.2]. One just needs to replace the gradient 
by the symmetric gradient. □ 

3.3. Decay estimate. In this section we derive decay estimates for our approximation h. 
The main ingredient is the following theorem which can be found in [11, Theorem 3.6]. It 
is valid in any dimension but needs <p" to be almost monotone. This is the only place in 
the paper, where we need this assumption on q>". 

Theorem 3.7. Let q>" be almost monotone. Ifh is a weak solution of (3.8), then there is 
an r > 2 such that for every ball Q <ZB with radius R > 

2 

R 2 ( f \VY(Dh)\ r dxj ' <cj\V{Dh)-(V(Dh)) Q \ 2 dx. 



The constants C and r depend only on the characteristics of <p and the constants in As- 
sumption 2.4. 

The regularity V £ W ' r with r > 2 ensures in two space dimensions that V is Holder 
continuous. This is the reason, why our estimates can only be applied to planar flows. It is 
an open question if V(Vu) is Holder continuous in higher dimensions. 

This provides the following decay estimates in the plane: 

Theorem 3.8. There exists y > such that for every X € (0, 1] 

j |V(Dh) - (V(Dh)) Xfl | 2 ^ < CX^j |V(Dh) - <V(Dh)) B | 2 ^-. 

XB B 

The constant C and y depend only on the characteristics of(p and the constants in Assump- 
tion 2.4. 

Proof. The result is clear if X > |,sowe can assume X G (0, A). Let R denote the radius 
of B. We compute by Poincare inequality on XB, Jensen's inequality with r > 2, enlarging 
the domain of integration and Theorem 3.7 

^|V(Dh) - (V(Dh)) AB | 2 c/x < C{XR) 2 j |VV(Dh)| 2 c/x 

XB XB 

2 2 

<C{XR) 2 (^j |VV(Dh)|'^ ' < CR 2 X 2 ^ (j |W(Dh)| r ^J ' 

xb i B 

<CA 2(1 -^^ |V(Dh) - (V(Dh)) B | 2 ^. 

B 

As r > 2 the proof is completed. □ 
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3.4. Proof of the main theorem. Theorem 1. 1 is a corollary of the following more general 
theorem. 

Theorem 3.9. Let B C M be a ball. Let u, it be a local weak solution of (1.1) on 2B, 
with <p and A satisfying Assumption 2.4. Let O) : (0,°°) — > (0, +°°) be non-decreasing such 
that for some /3 G (0, =r) f/ie function G)(r)r~P is almost decreasing, where y is defined in 
Theorem 3.8 and]? in (2.7). Then 

(3-9) |MI B MO B (B) + ll A ( Du )llBMO ffl (B) < cM L,2B( A ( Du )) + c ll G llBMO m (2B)- 

The constants depend only on the characteristics of <p and the constants in Assumption 2.4. 

Proof. The proof of the estimate of A(Du) follows line by line the proof of [12, Theo- 
rem 5.3]. It is based on Corollary 3.5, Lemma 3.6 and Theorem 3.8. We will therefore 
omit the proof here and restrict ourselves to the additional estimates for the pressure. 

We define H = A(Du) - G. It holds H G BMO^B) C BMO(B). We fix a ball Q C B. 
Then equation (1.1) implies that 

(3.10) V<i; G W '' 2 (£2) : {x- (7r) G ,div§) = (H-(H) e ,V§). 

Let E, G Wq' 2 (Q) be the solution of the auxiliary problem 

div<|; = 7t - (7t)q in Q, ^ = on dQ. 

The existence of such a solution is ensured by the Bogovskii operator [4] and we have 
||V^|| L 2(g) < C\\n — (^)g|| L 2(g). The constant C > is independent of Q. Inserting such t, 
into (3.10) we get 

II* - WqWIhq) = <*- We,div<^) = <H- (H) G , v^}. 

This and ||V§]| L2(e) < C\\n - (n) Q \\ L z {Q) implies \\n - (n) Q \\ L 2 {Q) < c ||H- (H) e || L2(G) . 
We find by Jensen's inequality 

(M^nf <j\K- (K) Q \ 2 dx <cj\H- {U) Q \ 2 dx < C||H|| BMO(Q) . 

e e 

In the last inequality we used the John-Nirenberg estimate, [15, Corollary 6.12]. It follows 
that % g BMO(B) and ||w|| BMO(Q) < C||H|| BMO(Q) . This implies that 

M ^Q {n) - c Wq) i|h|Ibm °(Q) - c H h IIbmo m (b) 

using the monotonicity of o. Since Q is arbitrary, we have IKHbmcWB) — I|H|Ibmo m (b)- 
Now H = A(Du) — G and the estimate for A(Du) conclude the proof. □ 

The choice (o(t) = 1 in Theorem 3.9 gives the BMO estimate. However, the choice 
0)(t) = t&, j3 G (0,27/7?') Theorem 3.9 gives the estimates in Campanato space J?? 1:2+ ^, 
compare [12, Corollary 5.5]. This is just Theorem 1.1. 

Remark 3.10. It is possible to transfer the Holder continuity o/A(Du) to Du and Vu. Let 

us discuss the case of power-law and Carreau type fluids. This follows from the fact that 
A -1 G for some o>0.IfK = 0, then a = rmn{l,p' - l}.IfK>0, then a = l. Now, 
A(Du) G C '^ implies Du G C '^ . Due to Korn's inequality we get Vu G C '^ as well. 
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Remark 3.11. Note that ifG G VMO(2B) in Theorem 3.9 we get that A(Du) G VMO(B). 
Indeed, since G G VMO(2B) there exists a nondecreasing function 6) : (0,°°) — > (0,°°) 
with lim r ^o &{ r ) ~ 0, such that ||G|| BMO (B r ) ^ ®( r )> f or a ^ Br C 2B. Defining (o(r) = 

a 

min{d)(r),rP T } we obtain by Theorem 3.9 the BMO m -estirnate for A(Du) and %, which 
imply that both are in VMO {compare [12, Corollary 5.4]). 

Remark 3.12. Let us now assume that the right hand side of (1.1) is not given in diver- 
gence form — divG with G symmetric, but rather as f G L s with s > 2. 

Let we W 2 - s {2B)(~}Wq S (2B) and aeW hs (2B) with {o) 2 b = be the unique solution of 
the Stokes problem — divDw + V<7 = f and divw = in 2B with w — on d (2B). Then f = 
-divG for G:=Dw- rrld and G is symmetric. Ifs = 2, then G eW l - 2 (2B) VMO(2B). 
Ifs > 2, then G G W Ls (2B) ^ «5f l - 2+( - 1 -^ (2B) = C ' 1 ^ (25). In particular, Theorem 3.9 
is applicable and for all s > 2 

N| _sf (j3) + HACDuJH^i^^ < ci?-^| B (A(Du)) + c||f|| LS(2B) 

/or s > 2 a«t/ J3 G (0, 1 — j] n (0, i^)- We additionally get VMO estimates if s = 2. 

The case s = 2 is obviously the limiting one in this setting. In the case of the p- 
Laplacian, i.e. no symmetric gradient and no pressure, it has been proven in [7, 16] that 
f G L" 1 (W 1 ) (Lorentz space; subspace ofL") implies A(Vu) G LT. It is an interesting open 
problem, if this also holds for the system with pressure and symmetric gradients (at least 
in the plane). Note that our results imply in this situation A(Du), 7t G VMO for n = 2. 

4. An application to the stationary Navier-Stokes problem 

In this section we present an application of the previous results to the generalized 
Navier-Stokes problem. We assume that u G W L f(Q.), divu = and K G if {SI) are local 
weak solutions of the generalized Navier-Stokes problem, in the sense that 

(4.1) V£ G W^{D.) : (A(Du),D«g> - (n,div^) = (G + u®u,D£) 

for a given mapping G : Q. > Rgym ■ 

In order to handle the convective term we need the condition 

(4.2) liminf > for some r > f . 

I-H-~ S r Z 

We have the following result 

Theorem 4.1. Let <p and A satisfy Assumption 2.4 and (4.2). Let u be a local weak solution 
of (4.1) on n. Let /3 G (0,=?) (a is defined in Theorem 3.8 and p in Lemma 2.7). IfBisa 

ball with 2B C ii and G G ^ l ' 2+ P(2B), then A(Du), n G jSf^+P (B). 

Proof. According to [1 1, Remark 5.3] we get that Du G L q (3B/2) for all q > 1. Conse- 
quently by the Korn inequality and the Sobolev embedding we get that u®u G Jz? 1 '" + ^(3B/2). 
Applying Theorem 1 . 1 we get the result. □ 

Exactly as in Remark 3.10 it is possible to transfer the Holder continuity of A(Du) to 
Du and Vu. 

Remark 4.2. A similar result has been proved also in [21], provided K > 0, by a completely 
different method, which requires the stronger assumption divG G L q (2B)for some q > 2. 

The same result was also proved in [32] for power law fluids with p G (3/2,2] and 
K>0, again under the stronger assumption divG G L q (2B)for some q > 2. 
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By our method we reprove these known results and improve them by weakening the 
assumption on the data of the problem. 

5. An Application to the parabolic Stokes problem 

Now we apply the previous results to the evolutionary variant of the problem (1.1). We 
set T > and / = (0, T), D. T = £i x / and assume that u G L°°(/,L 2 (£2)) with Du G LP(Q. T ) 
is a local weak solution of the problem 

d ( u-div(A(Du)) + V7T = f in£2r, 

(5.1) 

divu = in 

If the system of equations (5.1) is complemented by a suitable boundary and initial 
condition and if the data of the problem are sufficiently smooth it is possible to show 
existence of a solution that moreover satisfies 

(5.2) <?,uer°(/,L 2 (n)), 

see for example [23, 19, 6]. If we know such regularity of <9,u and f is smooth, it is easy to 
reconstruct the pressure % in such a way that n G L q (Q.T) with some q > 1 and 

(5.3) V^eCo(i2 r ): C + (A(Du) - Ttl, V£ )dt = 

Jo Jo 

The constant q is determined by the requirement A (Du) € L q (£lr). 

Applying the results from the previous sections of this article we obtain the next simple 
corollary. 

Corollary 5.1. Let A and (p satisfy Assumption 2.4. Let u G L C °(I 7 L 2 (Q.)) with Du G 
L^{Q.t) and divu = in £It solve the problem (5.1) and satisfy (5.2). Let B be a ball with 
2BcQ.andfe L°°(/,L 2 (n)). Then A(Du),7T 6 L°°(/,VMO(5)). 

Proof. The result is immediate consequence of <9,u G L C °(I 7 L 2 (Q.)) and Remark 3.12. □ 

Remark 5.2. Certainly, we can obtain a similar result for the problem (5.1) with con- 
vection, as soon as u (g)U G L°°(/,VMO(£2)). This follows for example from the fact that 
V(Du) G W i ' 2 (I,L 2 (Q,)) nL 2 (/,W 1,2 (£l)). Such kind of regularity is obtained, if it is pos- 
sible to test with d 2 u and Au. 

In [23] a method was developed to construct regular solutions of (5.1). The essential 
assumption was that the growth of A is sufficiently fast. It was necessary to assume that 

<p(s) A 

(5.4) liminf > for some r > |. 

i-H-oo S r J 

This assumption was not due to the presence of the convective term in the analysis of [23]. 
It was necessary to overcome problems connected with the anisotropy of the evolution- 
ary problem (5.1). The previous corollary is a first step to improve these results. If it is 
possible to show <9,u G L°°(I,L s (£l)) for some s > 2. Then for f G L°°(I,L s (£i)), we find 
by Remark 3.12 that A(Du) G L°°(I,C°>P (CI)) for J3 G (0, 1 - f- ] n (0, This implies 
(locally) bounded gradients Vu. So far the results of this paper are of local nature. An ex- 
tension of this technique up to the boundary would imply globally bounded gradients Vu 
and we could reconstruct the result of [23] for the generalized Stokes problem without the 
restriction (5.4). 
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